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ABSTRACT 


NASA's long range plan for the study of solar-terrestrial relations 
includes a Solar Probe Mission in which a spacecraft is placed in an 
eccentric orbit with perihelion at four solar radii. In the present 
concept of the mission, this orbit is attained by a Jupiter gravity- 
assist maneuver. A reasonable mission plan would require that there 
be no more than a 1% chance that the peak fluxes and fluences encountered 
near the Sun exceed those of the Jupiter flyby. It is shown that, if 
the assumptions in the calculations are correct, this requirement can be 
met. Some of the crucial assumptions underlying this conclusion are 
that, at four solar radii, there are no closed, stable configurations in 
which energetic particles are stored, and that energetic particles do 
not suffer appreciable energy loss in traveling from 4 solar radii to 
1 AU. The spacecraft could suffer fatal radiation damage (1) if it 
encountered one of the extremely large solar flares which occur, on 
the average, about once every ten years, (2) if there are invisible 
regions of trapped radiation and the spacecraft happened to pass through 
one of them, or (3) if there is some mechanism which causes an adiabatic 
energy loss so that the near-solar particles are appreciably more 
energetic than those observed at 1 AU. All these catastrophes are 
considered to be unlikely. It is also concluded that, although it is 
highly improbable that solar neutrons would be hazardous to the health 
of the spacecraft, there is a severe problem in limiting the neutron 
flux from a radioactive power supply to a level which allows solar 
neutrons to be detected. 
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I. 


INTRODUCTION 


NASA's long range plan for the study of solar-terrestrial relations 
includes a Solar Probe Mission in which a spacecraft is sent as close to 
the Sun as possible. In the present mission concept, the spacecraft 
components can probably be kept below their maximum operating temperature 
if the probe does not approach closer than 4 solar radii* (R^) to the 
center of the Sun. Figure 1 shows the near-solar portion of the planned 
traj ectory. 

The question arises whether or not the energetic particle environment 
encountered during such a solar flyby with perihelion at 4 might 
endanger the spacecraft and the mission. The Solar Probe Environment 
Workshop was formed to consider this problem. The Workshop met at the 
Jet Propulsion Laboratory on January 19 and 20, 1978. In addition to the 
authors of this report, the participants were R. W. Davies, T. N. Divine, 

L. A. Frank, E. F. Koprowski, J. E. Randolph, and E. C. Stone. 

Three different methods of placing limits on or estimating the 
energetic particle flux and fluence during the solar flyby are presented. 
First, in Section II, photon fluxes observed at 1 AU are used to derive 
upper limits for the proton and electron fluxes and fluences near the 
Sun. In Section III, we first use particle observations near 1 AU to 
estimate the maximum probability that the solar probe will encounter any 
given level of proton flux and fluence; maximum probabilities are ensured 
by making a wildly conservative assumption about the distribution of 
event sizes. We then make best guesses about the relation between the 
flux observed at 1 AU and the total number of particles released from 
the Sun and the distribution of event sizes to calculate a best-guess 
probability that the solar probe will encounter any given level of 
proton flux and fluence. 

Electron and neutron bombardment are discussed in Sections IV and 

V, respectively. The overall outlook for a mission with perihelion at 

4 R is discussed in the final Section, 
s 


* 1 R g £ 0.005 AU £ 7 x 10 5 km 
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II. UPPER LIMITS ON ENERGETIC PARTICLE FLUXES AT THE SUN FROM 

PHOTON OBSERVATIONS 

Energetic electrons and ions produce bremsstrahlung x-ray emission 
and nuclear gamma-ray line emission, respectively, through collisions 
with the solar atmosphere. The electrons also interact with the ambient 
plasma and magnetic field to produce radio emission. Although the 
mechanisms for coronal radio emission are not quantitatively well under- 
stood, ground-based and satellite observations can provide information 
on the location and extent of the energetic particle populations to 
distances of many solar radii (Fainberg and Stone, 1974). Observations 
of the x- and gamma-ray emission can provide quantitative estimates of 
the spectrum and fluxes of energetic particles in the emission region. 
These estimates can be used to obtain essentially model- independent 
upper limits to the energetic proton and electron fluences for the Solar 
Probe Mission. 


Following Lin (1974), consider a volume, V, containing ambient 
plasma and energetic particles. Let 


E 

hv 

n H 

dn(E)/dE 

da(E,hv)/d(hv) 


= kinetic energy of an energetic particle with velocity v 
= energy of photon 

= number density of the hydrogen nuclei 

-3 -1 

= density of energetic particles (cm keV ) 

= differential cross section for photon production 

f 2 . - 1 , - 1 . 

(cm xon keV ) . 


Then the photon flux at 1 AU, assuming the particles are isotropic, is 


dJ (hv) = n H V f do dn 
d(hv) ’ 4*R 2 i V d(hv)dE 


( 1 ) 


13 

where R = 1 AU = 1.5 x 10 cm. For non-relativistic electrons the Bethe- 
Heitler cross-section is applicable 


da(E,hv) 

d(hv) 



( 2 ) 


Equation (1) can be inverted with (2) to obtain the electron spectrum. 
For a power law x-ray spectrum 
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dJ/d(hv) = A(hv) ^ 


(3) 


we obtain 


dJ t 

dl" 


vdn 

dE 


2.27 x 10 51 y (y-1) 2 B (y-%, 3/2) AE 


-y+1 


n H v 


(4) 


where B(x,y) is the beta function. 

The strongest solar y-ray line, both from theory (Ramaty et al., 
1975) and observation (Chupp et al., 1973), is located at 2.22 MeV. This 
line results from neutron capture by hydrogen and the subsequent 
deexcitation of the deuterium. The neutrons are produced by energetic 
protons and heavier ions interacting with various nuclei. Detailed 
calculations have been performed on the neutron production and propagation 
to obtain the 2.22 MeV photon yield (Wang and Ramaty, 1974). For typical 
power law proton spectra with exponents of -2 to -3, the 2.22 MeV flux at 
1 AU is given by 


J 2 2 (photons cm 2 sec "S ^ 5 x 10 ^n (>io MeV)n fI V 


(5) 


Thus 


Jp(>10 MeV) (cm^sec) ^ 


■v 


n v 
P 


'U 

'V 


10 


55 


u H V 2.2 


( 6 ) 


It is important to note that these relationships hold at a given 
instant of time. That is, the instantaneous number and spectrum of 
energetic protons or electrons in the region is related to the 
instantaneous gamma-ray or x-ray emission. 

Most of the x-ray and y-ray emission probably originates from 
regions much closer to the Sun than the solar probe will traverse. An 
absolute upper limit to the energetic particle fluxes would thus be 
obtained under the assumption that all the emission comes from the 4 to 
10 R region. The density n^ estimated from radio observations (Fainberg 

® . U r 

and Stone, 1974) is ^3 x lO^cm - ^ at 10 R and ^3 x 10 cm at 4 R . We 
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assume the particles fill a region steradian wide and ^2 R thick at 

34 3 s 

4 R (total volume ^10 cm ) and that the volume at other distances 
s 3 OQ 

scales as r . Then n V ^ 4 x 10 J , approximately constant from 4 to 10 R . 

H s 

Taking the observed x-ray and y-ray emission (or upper limits) for quiet, 

one-a-day flares, and very large flares, and using equation (4) and (6), 

we have computed the maximum particle fluxes at 4 to 10 R distances 

s 

(Table 1) . The total integrated particle fluences depend on the flare 

2 3 

emission period, At £ 10 -10 sec. For non-flare fluxes. At £ 10 sec, the 
traversal time across 1 steradian solid angle. 

These flux and fluence limits are quite high, but more sensitive 
observations with better x-ray and y-ray instruments over the next few 
years can substantially reduce these limits. 


TABLE 1. UPPER LIMITS ON ENERGETIC PROTONS AND ELECTRONS FOR SOLAR PROBE 



Quiet-time 

Typical 

Very large 

I. Energetic protons 

(non-flare) 

flare 

flare 

2.22 MeV photon flux* 

<5x10" 3 

<5xl0~ 3 

'vO. 3 

, 2 

(cm sec) 


<10 13 

<io 15 

>10 MeV proton flux 

<10 13 

, 2 x~l 

(cm sec) 

<io 17 

<io 16 

<io 18 

>10 MeV proton fluence 

-2 

cm 

(At = 10 4 sec) 

(At = 10 3 sec) 

(At = 10 3 sec) 


II. 


Energetic electrons 
X-ray parameters** 


dJ 

d (hv) 


= A(hv)" Y 


A 

Y 


>20 keV electron flux 

, 2 v"l 

(cm sec) 

>1 MeV electron flux 

, 2 x"l 

(cm sec) 

>20 keV electron 

-2 

fluence (cm ) 

>1 MeV electron 

-2 

fluence (cm ) 


<20 

'VjxlO 3 

VLxlO 6 

^2 

^4.5 

n,2.7 

13 

<2x10 

<3xl0 14 

<2xl0 17 

— 

<2xl0 10 

<io 14 

<2xl0 17 

<3xl0 16 

<2xl0 20 

(At = 10 4 sec) 

(At = 10 2 sec) 

(At = 10 



12 

<2x10 

<io 17 


*From Chupp et al. (1968) and Chupp et al. (1973). 

**From Datlowe et al. (1974), Lin and Hudson (1976), and Kane (1978). 
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III. PROTON FLUXES ESTIMATED FROM PROTON OBSERVATIONS 

A. ASSUMPTIONS 

Protons and electrons are the most abundant solar energetic 
particles detected in interplanetary space and are the most likely to 
be observed at a few solar radii. 

At 1 AU, the observed solar particle peak intensity ranges from 
3 2 

greater than 10 particles/cm -sec-ster at energies above 20 MeV down 
to micro events with integral intensity above 20 MeV some five orders 
of magnitude lower. The energies of particles range from tens of keV 
up to some 20 GeV in the infrequent high energy solar particle events. 

Below 20 MeV, a variety of different types of proton events has 

been observed, such as flare associated particle events, corotating 

particle streams and energetic storm particle events in addition to a 

small galactic component. The origin of many of these event types is 

now becoming well understood. The flare associated events have a 

well defined solar signature starting with an active region, radio 

and x-ray emission, and a solar flare. In the case of the corotating 

streams of energetic particles which represent the dominant component 

observed during solar minimum, there is strong evidence of a positive 

gradient between 0.3 AU and 3-5 AU (McDonald et al. , 1976; Shea and 

Smart, 1973; Simnett, 1972; Van Hollebeke et al., 1975; and Van Hollebeke 

et al., 1978) which indicates that acceleration is taking place in the 

interplanetary medium. For those events, it is therefore expected that 

the intensity of low energy particles will be smaller near the Sun than 

2 

in interplanetary space (typically £ 10 particle/cm -sec-sr-MeV at 
1 MeV at 1 AU with an exponential momentum spectrum of the form 
dJ/dP « exp(-P/P Q ) with P q = 9-15 MV and thus of small importance 
compared to the intensity of flare associated particle events. 

Since the observations of protons below 20 MeV are greatly 
influenced by interplanetary propagation effects (acceleration, as 
indicated above, deceleration, diffusion, shock effect, etc.), it is 
not easy to infer the solar condition at those energies from inter- 
planetary observations. 

However, at higher energy, these interplanetary effects are less 
important, thus making possible the use of proton observations above 
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20 MeV to define energetic particle conditions at the Sun. Note also 

2 

that since the range of a 20 MeV proton is ^0.5 g/cm or 0.2 cm in 
aluminum (it is ten times smaller at 5 MeV), this energy range is also 
the most relevant one in defining the expected radiation environment 
which could be hazardous to a spacecraft with perihelion between 2 and 
10 solar radii. However, if needed, the conclusion drawn from proton 
observations above 20 MeV can be extrapolated to low energy. 

In estimating the probabilities of encountering different levels 
of radiation, we make three basic assumptions: 

i) We assume that there are no permanent storage configurations at 

altitudes as great as 4 solar radii. The evidence that stable closed 

loops in the corona do not extend beyond 2.5 R g is rather extensive. 

Observations of the eclipsed Sun allow the density structure of the 

corona to be established out to about 4 or 5 R g . Assuming the density 

structure to outline the coarse configuration of the field, we note 

that reviews (Newkirk, 1967) of such observations covering a century of 

eclipses set the upper limit of stable closed coronal forms at 

R 'u 2.5 R . Infrared and satellite observations out to 13R/->, confirm 

this conclusion (Tousey et al., 1974). Theoretical calculations of 

the interaction of coronal magnetic fields and the expanding solar wind 

plasma have grown quite sophisticated in the last decade (see Pneuman, 

1968 and 1973) and reach the same conclusion. Above ^2.5 R the dynamic 

s 

forces of the solar wind are sufficient to open up the magnetic field 
into interplanetary space. 

Some investigators argue that invisible, closed, stable magnetic 
loops (ICSML) may still exist above the canonical 2.5 R g . Although 
caution must always be exercised in confronting the invisible, the 
proponents of such ICSML’ s have thus far failed to answer: 

• How do ICSML’ s avoid containing material and thus evade detection 
while other quite substantial coronal structures are well documented 
in the region 2.5 - 13 R g ? 

• How do ICSML’ s avoid being swept out by the solar wind? 

• What phenomena indicate their existence? 
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ii) We assume that at a given radial distance above 4 R g the total number 
of particles at a given energy, crossing unit area during a sufficiently 
long time, is the same at all solar latitudes and longitudes. Flare 
events occur at many different locations on the Sun, and energetic 
particles can propagate through the corona many tens of degrees in 
longitude and also presumably in latitude. At least within solar-active 
zones( e.g., <30° in latitude) the total number of emitted particles should 
thus be spread evenly in longitude and latitude. However, this assumption 
may cause us to overestimate the fluences at high solar latitudes, which 
are far from active regions. 

iii) We assume that energetic particles propagating from 4 R g to 1 AU do 
not suffer appreciable energy loss; i.e., a 20 MeV proton observed at 
Earth was a 20 MeV proton at 4 R^. Any energy loss would result principally 
from adiabatic deceleration in the expanding solar wind. However, particle 
mean-free paths in the interplanetary medium appear to be too large for 
this to be a major effect, except perhaps during the decay with time, at 
late times during flare events. 

B. MAXIMUM PROBABLE PROTON FLUX AND FLUENCE 

According to Reedy (1977), in the eight years between 1965 and 1973 

the total fluence at Earth of solar protons with energy >_ 20 MeV was 

9 -2 -2 -1 

3 x 10 cm , which corresponds to an average flux I = 12 cm s 

avg 

If the flux were organized into discrete events, each with intensity I, 

then an event would be in progress a fraction I /I of the time. Since 

avg 

the maximum probability of observing a given flux level I would correspond 
to the assumption that all of the solar events during this 8 year period 
had this same intensity I, we find that 


P (I) at 1 AU = I /I. 
max avg 


If we further assume that the solid angle of the event is independent of 
distance from the Sun, then P max (^) would vary as the inverse square of 
the solar distance, and, at 4 solar radii, 


P (I)--P 
max I 


4 R 
s 

]. AU 


3 x 10 


( 7 ) 


-2 -1 

for protons above 20 MeV, where I is in units of protons cm s 
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Clearly, (7) is a simple result. It yields, for example, that 

there is a maximum of a 1% chance that in a given region near 4 R the 

6 2 s 

spacecraft will encounter a flux of 3 x 10 protons/cm -sec, with 

energies >_20 MeV. Since this result holds for any region near 4 R g 

through which the spacecraft passes, there is also a maximum probability 

of 1% that it will encounter a fluence of It , where t is the time 

sp ^ sp 

spent by the spacecraft near perihelion. With t ^16 hours, this 

11 2 S P 

fluence is ^2 x 10 protons/cm . 

C. MOST LIKELY PROTON FLUX AND FLUENCE 

We next attempt to calculate a more realistic probability of 
encountering given proton fluxes and fluences by inferring the total 
number and energy spectrum of particles released at the Sun in each 
event as well as the probable frequency of events of different sizes. 
Because of the many assumptions and extrapolations made, the resulting 
probabilities may be less reliable than the upper limits and maximum 
probabilities calculated in the two previous Sections. 

1. TOTAL NUMBER OF PARTICLES INJECTED AS A FUNCTION OF 
OBSERVED PEAK FLUX 

There are several methods of estimating the total number of particles 
released into the interplanetary medium during a solar flare event. We 
first make the assumption that we are only interested in particles that 
reach VI AU (i.e., we assume that essentially all particles that 
propagate out to perihelion of the Solar Probe orbit will eventually 
reach 1 AU without significant energy loss). Then two methods of esti- 
mating the total particle injection are attractive. Spherically symmetric 
diffusion models with impulsive injection give a simple result for the 
total number of particles that depends only on the peak flux and the 
radius at which it is observed. However, this method is model dependent 
and the specific assumptions of impulsive injection and spherical 
symmetry are in direct contradiction to observations. 

Another technique for estimating the total particle injection can 
be constructed from the continuity equation. Essentially, the net 
energetic particle streaming along the flux tube of the spacecraft may 
be computed from the observed flux and anisotropy measurements. The 
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net streaming need just be summed throughout the event to determine the 
total number of particles escaping through the flux tube. This technique 
has the advantage that it is independent of any model of interplanetary 
energetic particle propagation and relies purely on observed quantities. 
However, there are difficulties in implementing this technique since it 
requires much more data than just the peak flux which is all that is 
needed by the spherically symmetric diffusive approximation with impul- 
sive injection. For the net streaming technique both the flux and the 
first harmonic of the Fourier decomposition of the anisotropy must be 
known as a function of time throughout the event with sufficiently fine 
time resolution and statistical accuracy to define the outward current. 
Unfortunately, there have been very few detectors with angular sectoring 
and sufficient geometric factors and duty cycle to make statistically 
significant anisotropy measurements in the 20 MeV range; existing 
measurements primarily cover only a few of the larger events during 
the latter part of solar cycle 20. However, high quality measurements of 
20 MeV proton time histories are available during most of ,the solar cycle. 
We show below that the particle streaming in the flux tube of the space- 
craft may be estimated to within about a factor of 2 purely from the flux 
time history thereby removing the major problem with the implementation 
of this technique. But first we shall apply both the diffusion model and 
the streaming technique to a sample solar particle event. 

a. APPLICATION OF THE BASIC TECHNIQUES 

We shall apply both the diffusion model and the streaming technique 
to the October 30, 1972, solar particle event. This event has been 
chosen because of the availability of anisotropy data from IMP-7. It is 
one of the largest events during solar cycle 20 and is also one in which 
the time histories at >10 MeV look "classical" and therefore the two 
methods might be expected to yield similar results. However, the time 
span of interplanetary particle injection at lower energies was quite 
extended during this event and if there was also an extended injection at 
>20 MeV, it could contribute significant disagreement between the two 
methods. 

The data are from the 25.2 - 49.5 MeV proton channel of the APL/JHU 
experiment on IMP-7 which has 8 angular sectors in the ecliptic plane. 
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Three hour averages of the differential flux at 37.3 MeV and the first 
harmonic of the angular distribution are listed in Table 2 for the 
period when the flux was greater than 10 (cm sec sr MeV) which is 
within about two orders of magnitude of the peak flux. Let us now 
look at the total particle release estimates for an observer at 1 AU. 

i. SPHERICALLY SYMMETRIC DIFFUSION MODEL WITH 
IMPULSIVE INJECTION 

The diffusion model gives the simple result 


N 


8J r 3 
m 

v 



-3/2 


( 8 ) 


where N is the total number of particles /MeV, J is the maximum 

m 

observed differential flux, r is the distance from the sun and v is the 
particle velocity. At 1 AU this reduces to 

N = J x 5. 3 x 10 30 . 
m 

Substituting the peak flux from Table 2 


29 -1 

N = 5.1 x 10 MeV . 


TABLE 2. FLUX AND ANISOTROPY DURING THE OCTOBER 30, 1972 EVENT 


25.2-49.5 MeV FLUX ANISOTROPY 


DAY 

HOUR 

(10-3) 

MAGNITUDE 

304 

9-11 

6.1 

.92 


12-14 

11 

.26 


15-17 

17 

.38 


18-20 

28 

.21 


21-23 

48 

.17 

305 

0-2 

85 

.36 


3-5 

57 

.34 


6-8 

28 

.16 


9-11 

16 

.17 


12-14 

11 

.13 


15-17 

5.4 

.25 


18-20 

1.6 

.42 

Peak flux 9.7 

-2 -2 - 
x 10 cm s 

’ 1 sr _1 MeV -1 at 0100-0200 

on October 31. 
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ii. OBSERVED STREAMING 

If energy changes are insignificant, the total number of particles 
flowing out along a flux tube during a flare event is 



where S is the net particle flux and A is the area of the flux tube. 

Early in the event, when the anisotropies are large and field aligned, 
the volume being considered is truly a flux tube and S ^ Sy. The net 
current is essentially the measured anisotropy as long as the magnetic 
field vector is close to the ecliptic plane. Starting near the peak of 
the event, when the anisotropies usually cease to be strongly field 
aligned, the radial component of the observed anisotropy may be taken as 
a good approximation to the particle loss. 

Let us now apply this technique to the October 1972 event assuming 
spherical symmetry so that the result may easily be compared with the 
diffusive model calculation above. Then 

N = 4tt (1 AU) 2 ZJ.S.T. (9) 

where is the observed flux during the ith time period t , and is 
the effective outward streaming anisotropy. Using the observed 3-hr 
anisotropy amplitudes for we get 

30 -1 

N = 2.7 x 10 MeV . 

This is a factor of “v5 greater than the estimated total particle release 
from the spherically-symmetric, impulsive- inj ection diffusion model. 

Since this calculation has also assumed spherical symmetry, it's the 
impulsive injection assumption in the diffusive model that has probably 
contributed the greatest source of error. 
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b. ANGULAR DISTRIBUTIONS 

Clearly spherical symmetry is not what is observed; however, there 
have been very few observations of the angular dependence of solar 
particle events. The cleanest observation reported has been the 
April 10, 1969, event which was observed by 5 spacecraft covering 
more than 180° of heliographic longitude (Gold et al. , 1977a). The 
intensity of > 10 MeV protons during this event was found to fall off 
exponentially from the center of a preferred release zone with an 
e-folding angle for these azimuthal gradients of v25°. The center of 
the distribution during this event, however, was ^100° west of the 
flare site. Statistical studies of the peak flux during an event as a 
function of the angular distance to the flare site have been presented 
by McKibben (1972), Reinhard (1975), and Van Hollebeke et al. (1975). 

They confirm the existence of azimuthal gradients. However corotation, 
possible displacements of the distribution center from the flare site, 
etc., preclude accurate estimates of the gradients from statistical 
studies. If we use the distribution of peak fluxes for events associated 
with flares from 20-60° west and 20-60° east solar longitude, as reported 
by Van Hollebeke et al., we derive an estimate of e-folding angles of 
25 to 70°. McKibben's observations of the May 1968 event by IMP-4, 
Pioneer 6, and Pioneer 7 lead to e-folding angles of 30-50° at M.5 MeV. 
Gold et al. (1977b) and Roelof et al. (1975) have examined the azimuthal 
gradients in individual events by solar wind mapping techniques. They 
found e-folding angles ranging from 10 to 30° with an average near 20° 
that were approximately energy independent; however their work was in 
the 0.3-5 MeV range and the gradients may not have the same values at 
£ 20 MeV. 

In spite of the limited direct evidence, we shall assume that all 
solar particle events in the £ 20 MeV range have azimuthal gradients 
with e-folding lengths of ^20-30°. Then, normalizing to the center of 
the distribution, the ratio of the total number of particles injected 
into the interplanetary medium in a distribution with e-folding angle 
9 q to the total injection in the spherically symmetric distribution is 



-13- 


N(0 ) 1+ exp (-n/ 6 ) 

o _ o_ 

N (spherical) 2 (1 + 1/0 3 ) 

o 

where 0 is the e-folding angle. Thus, with gradients, the total injections 
o 

are only ^0.05 to 0.1 of the spherically symmetric estimates. For the 
example above, applying the observed streaming technique to the October 1972 
event, an azimuthal gradient with an e-folding length of 30° would, if the 
Earth were at the center of the distribution, reduce the estimate of the 
total injection to 


N 


29 -1 

x 10 MeV 


Since we are not certain that the Earth was exactly at the center 

29 -1 

of the distribution, we will assume that N = 4 x 10 MeV for the 
October 1972 event, for which the peak flux was 10 ^ cm ^sr ‘'"MeV ^ . 

If we follow Equations (8) and (9) , assume that the total number of 
particles injected is proportional to the peak flux, and use the October 
1972 event for normalization, we obtain 

N = 4 x 10 30 J (10) 

m 

for J in units of cm 3 s ^sr ^MeV ^ and N in units of MeV 
ra 

c. PROTON ENERGY SPECTRA 

The value of N calculated above is the total number of particles 
injected per MeV. To calculate N q , the total number of particles 
above 20 MeV, we must make some assumptions about the proton energy 
spectrum. 

The energy spectra observed at 1 AU represent the combined effects 
of the source spectrum, coronal transport, the release into the interplanetary 
medium, and interplanetary diffusion and velocity dispersion. It has 
been shown that in measuring the particle spectrum near the maximum 
particle intensity, interplanetary effects are minimized. This spectrum 
is then representative of the energy spectral distribution both close to 
the Sun and after release into the interplanetary medium (Van Hollebeke et al., 
1975). 
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Based on observations of some 100 solar flare-associated particle 
events detected at 1 AU, it was found that for protons in the energy 
range 20-80 MeV, a power law in differential kinetic energy of the 
form J « E ^ appears to give the best representation. Figure 2 
shows a representative sample of spectra. The spectral index y varies 
with the flare longitude relative to the observer from y = 2 to y = 5.5 
as shown in Figure 3. It has been argued that for heliolongitude 20-80° 
west of the Earth's central meridian (preferred connection region for an 
observer near Earth) the spectral index y is representative of the 
source spectrum. For these heliolongitudes, y displays a surprisingly 
small range, varying between 2.0 and 3.2 and a mean of 2.7. The 
steepening of the spectra with increasing angular separation from the 
preferred connection region can be explained in terms of an energy or 
velocity dependent particle escape rate from the corona. 

It is therefore expected that at any given point outside the 
escape boundary into the interplanetary medium, the spectral distri- 
bution of the particles will be dependent on the observer's position 
relative to the flare site: near the flare site, the particle population 

will display the characteristic source spectrum: 


- 2.1 + 0.5 


(ID 


while away from the flare, the spectrum will be poor in energetic 
particles as those will have escaped faster (thus, a steep spectrum). 
From the point of view of modeling the expected particle environmental 
condition at a few solar radii, it is reasonable to assume an average 
spectrum equal to that of the characteristic source spectrum. 


N, the total injection per MeV, was determined above from measure- 
ments over the energy range 25.2 to 49.5 MeV. Equation (11) can be used 

to convert N to N as follows : 
o 


N 

o 


(49.5 - 25.2) N /e 
20 



25.2 


dE 


53 N 


or, 


32 

N = 2 x 10 J . 

o m 


( 12 ) 
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2. FREQUENCY OF EVENTS VERSUS PEAK FLUX 

a. ENERGETIC SOLAR PARTICLE EVENT PROBABILITY 

The frequency of particle events produced at the Sun can be readily 

deduced from the nearly continuous interplanetary observations made with 

Earth orbiting satellites at 1 AU during solar cycle 20 (1964-1976). 

Figure 4 shows the number of > 20 MeV proton events with peak intensity 
2 

>0.5 protons/cm -sec-ster detected at 1 AU between 1962 and 1975. Also 
shown in the figure is the number of sea level neutron monitor events 
(with energy > 500 MeV/nucleon) from 1956 to 1972 (Shea and Smart, 1973; 
McDonald, Fichtel and Fisk, 1974). The solar cycle dependence of event 
frequency can be easily appreciated for > 20 MeV proton events, using 
the Zurich sunspot number as a simple index to indicate the overall 
level of solar activity. For solar cycle 20 (which represents a good 
approximation of the mean solar activity between solar cycle 8 and 20, 
as indicated in Solar Geophysical Data, 1977) there is a linear relation 
between the number of events observed per year and the yearly- averaged 
sunspot number as shown in Figure 5. In this figure, the error bars 
attached to the number of events represent a lo statistical fluctuation. 

Note that the number of events/year represents all the > 20 MeV 

2 

proton events with peak intensity >0.5 protons/cm -sec-ster detected 
by Earth orbiting satellites during every year from 1964 through 1975. 
Assuming that the event visibility conditions are the same at a few 
solar radii as at 1 AU, the event frequency seen by a solar probe during 
perihelion passage should be the same as that indicated in Figure 5. 

The event visibility at 1 AU has been interpreted in terms of coronal 
transport and release into interplanetary space. It is therefore 
expected that the event frequency outside the corona (> 3 R g , below 
which coronal transport may take place) will be the same as at 1 AU. 

The equation of the straight line through the data in Figure 5 
is: 

Number of events per year = 0.15 R - 0.60 

where R^ is sunspot number. The probability of an event occurring in 
any hour is thus given by 

P(>0.5 cm _2 s _1 sr _1 MeV _1 ) = (1.7 R z - 6.8) x 10~ 5 . 


(13) 
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b. SIZE DISTRIBUTION OF FLARE ASSOCIATED PARTICLE EVENTS 

The energetic solar particle event frequency is not only a function 

of time but it is also a function of intensity and energy. The number 

of events per unit intensity as a function of intensity is displayed in 

Figure 6 for all the events detected by the GSFC cosmic-ray experiments 

on the Earth orbiting satellites IMP 4 and IMP 5 between 1967 and 1972 

(Van Hollebeke et al., 1975). To minimize the effect of the spectral 

dependence upon longitude, the peak proton intensity measured at 40 MeV 

is used to characterize the size of the event. The galactic background 

observed at 1 AU fixes the lower limit of detectability of event 

intensity to vLO ^ particles cm ‘'"sr "'‘MeV ^ above 20 MeV. The size 

distribution of. all detected particles can be adequately described by 
“Ct 

dn/dJ « J where J is again the maximum differential intensity at 

m m - 2-1 -1 -1 
1 AU, in units of particle cm s ster MeV and dn is the number of 

events in the range of intensity dJ. a is found to be 1.15 + 0.05. 

The size distribution of the flare- identified events is also plotted 

in this figure. Note that the form of the distribution is the same for 

the flare-identified events and for the total number of events. It is 

also found that the form of the distribution is independent of the 

observer's position relative to the associated flare except that the 

number of events per unit of intensity decreases by a factor n<2 when 

the observer's position relative to the associated flare is > 80° away 

from the best connection region. Note that in Figure 5, the number of 

2 

events per year for > 20 MeV proton with peak intensity > 0.5 protons/cm s 

-3 2 

is the fraction per year of the integral above ^6 x 10 protons/cm -sec- 

- 2. 7 

sr-MeV obtained from Figure 6 assuming J « E . 

When we combine this result with Equation (13) we find that the 
probability of observing a peak flux >_ J q in any one hour is given by 

30 I oo ^ ~ 0 • 15 

J 6x10“ 3 

o 


P(> J Q ) = 


(1.7 R - 6.8) x 10 


-5 


/ 


sr 


The hourly probability of the injection of > N particles can then be 

— o 

found from Equation (12) : 
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P (>N ) 
— o 


(0.55 R z 


2.2) N 

o 


-0.15 


(14) 


3. EXPECTED PROTON FLUX AND FLUENCE 

We are concerned with the encounter of the Solar Probe with 
energetic particles, particularly as they influence the two parameters: 


I 



J(0,(j)) dm dE 


0 E 0 


the mean intensity above some cut-off 
energy (E^ 20 MeV for protons) 


and 


F 



the fluence of such particles 


During its perihelion passage at 4 R the Probe will spend ^ 16 hrs 

4 s 

(5.8 x 10 sec) in the region R < 10 R and may encounter energetic 

s 

particles of several distinct populations: 

(1) Particles escaping on open field lines on their way into 
interplanetary space, 

(2) Particles stored in stable closed magnetic configurations, and 

(3) Particles temporarily resident in closed, rapidly evolving 
magnetic configurations. 

The distinction between the three cases will be examined below. 
Throughout we shall attempt to define the probability that a specific 
event will befall the Probe during its 16 hours of passage below 10 R g . 
The following illustration is useful. 

From the last subsection, we know P(N q ), the probability that in 
any given hour an event will occur with a total particle output N q 
above ^ 20 MeV for protons. The probability that during 16 hours at_ 
least one such event will occur is 
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P'(N ) = 1 - (1 -tP(N )] t ) where t = 16 

o O * 


and for P(N ) << 1 
o 


P'(N ) = 1 - 1 + tP(N ) + . . . £ tP(N ) 

o o o 

a. CASE (a) (IMMEDIATE ESCAPE) 


At a distance R we have the approximation 


I 

a 


N 


2 

UjR T 


where id = solid angle of the escape region as seen from the 
center of the Sun 


t = event duration. 


Also 


F 'v I — 
a a v 


sc 


where 

l = length of spacecraft trajectory in the region w 
v gc = spacecraft speed. 

The probability of the spacecraft’s encountering such an event in any 
hour is 

P(I a ) 2 P(N 0 ) 


and that at least one such event be encountered during perihelion passage 
is 
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If we take w = 1 sr, R=8R =5.6x 10^ cm, t = 16 hr = 5.8 x 10^ s, 
7 -1 S 

and v = 3 x 10 cm s , then 
sc 


-29 -2 -1 

I =5.5x10 Ncm s 
a o 


( 15 ) 


-24 -2 

F = 1.0 x 10 N cm 

a o 


(16) 


P' = 77 P(N ) 
4tt o 


1.27 (0.55 R_ - 2.2) N 
Z o 


-0.15 


The probabilities P' of encountering flux levels I and fluences F due 

a a 

to immediate escape for years in which R = 100 and 200 are plotted in 

Z 

Figures 7 and 8 together with the maximum probabilities calculated in 
Section IIIB. 

b. CASE (b) (STABLE CLOSED LOOPS) 

As discussed in Section IIIA, we believe that there are no stable 

closed loops which extend to 4 R . Therefore P'(I, ) = 0. 

s b 

c. CASE (c) (CORONAL TRANSIENTS) 


For this case we assume that all particles which escape are 
itaneously inject 
particle distribution 


instantaneously injected into a volume ^ trans - Then for an isotropic 


. N v 

T „ 1 O 

I “V — 

c 3 V 

trans 


and 


F 'v I — 

C - C V 


sc 


where v is the particle speed. For the probability, we take 


P(I ) = P T 0 (N ) f- 

c T o 4 tt 


( 17 ) 


where P is the probability of a transient event anywhere on the Sun 
during any hour, 0(N q ) is the fraction of transients with particle 
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output > N , and oj is the solid angle subtended by an approximately radial 
— o 

transient. To estimate P_ 0(N ), we can either assume (i) that the only 

T o 

transients which contain energetic protons (>_ 20 MeV) are those for 
which energetic particles are detected at 1 AU, in which case 


P_0 (N ) = P (N ) (18) 

To o 

where P(N ) is given by Equation (14), or (ii) that most transients 
o 

contain energetic particles which do not reach 1 AU. Since transients 
have not been observed to return to the Sun (Gosling et al. , 1974) and 
a large fraction of particles injected into interplanetary space should 
eventually be detected, we conclude that the second case is unrealistic 
and that only a small fraction of transients contain energetic particles. 

We note, however, that since the average rate of coronal transients 
during 1973 was 0.5 day ^ (Hildner et al., 1976), whereas Figure 4 shows 
6 energetic proton events in 1973, only 3% of the transients led to a 
detectable proton event. If all transients do contain energetic 
particles and the distribution of event sizes is the same as at 1 AU, 
then we may be underestimating the radiation hazard by a factor of 30; 
we believe this to be unlikely, however. 

From Equations (17) and (18), the probability of at least one 
encounter with a transient containing energetic particles is 


P' (I ) = t(hrs) P(N ) y— 

c o 4 it 


For the geometry of the transient we consider two models - that 

of a "bubble" and that of a single flux "rope." We assume that a bubble 

occupies a cone of 60° as seen from the center of the Sun and extends 

from 4 to 10 R . Then w = 0.84 sr and V = 245 R ^ = 8.3 x 10^ cm^. 

s 7 trans 1Q s_ 1 

With 2. = 8 R , v =3x10 cms , and v = 10 cm s , we obtain 
s sc 


, n , n -26 -2 -1 

4.0 x 10 N cm sec 

o 

-22 -2 

7.5 x 10 N cm 

o 


p’d c ) 


- 1.1 (0.55 R_ - 2.2) N 

Z o 


-0.15 
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For a rope transient of height h and rope radius a, 


and for 


Then 


V 'v 3h Tra^ 

trans 


h = 10 R 

s 


a 

= 0.2 R 

V 

s 

33 3 

= 1.27 x 10 cm 

trans 

£ 

= 2a. 


o * in" 24 xt “ 2 "I 

2.6 x 10 N cm sec 

o 

-21 -2 
2.3 x 10 N cm 

o 


oo 

4tt 


2ra 


4tt ( 8 R )‘ 
s 


= 3.1 x 10 


-4 


P'(I ) = 5 x 10 -3 (0.55 R - 2.2) N ° -15 
c L o 


The probability of encountering any flux or fluence levels due to 

the Probe passing through a bubble-shaped transient is also plotted in 

Figures 7 and 8 for a year in which the sunspot number is 200. The 

ratio of the probability of encountering a given flux or fluence due 

to an immediate escape event to the probability of encountering the 

same flux or fluence due to a bubble-shaped coronal transient is 

P' (I )/P'(I ) = 1/3. The probability of encountering a rope shaped 
d c 

coronal transient is negligibly small (it is below the bottoms of 
Figures 7 and 8). 
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IV. ELECTRONS 

A. OBSERVATIONS 

To a first approximation, the distribution of the number of 
electrons at the Sun can be deduced from the distribution of protons 
calculated from the parameters described above, at least for high 
intensity events, which are the main concern of this study. This 
approximation is based on two results: 


la) 


From Equation (11) , it follows that the characteristic 
source density spectrum (i.e., the number of particles per 
unit energy and volume) for protons between 20-80 MeV is of 
the form 


6 (E) « 
P 


(!) 


-3.2 + 0.5 


and may be assumed to be of the same form down to 10 MeV. 
At even lower energy, the spectrum may be the same or 
flatter. 


lb) 


The characteristic density spectrum for electrons had been 
reported as 


6 (E) « 
e 



and 6 (E) « 
e 


-3 + 0.2 



between 3 and 12 MeV (Simnett, 1972) 
between 12 and 45 MeV (Datlowe, 1971). 


Below 3 MeV, a number of observed spectra seem to fit a 
power low in energy with spectral index 2 to 4 (Lin 1974) . 

2) The ratio between the 0. 2-1.1 MeV electron maximum intensity 
and the 10 MeV proton maximum intensity was found to be 
constant for 80% of the events observed at positions well 

-2 

connected to flare regions and with peak intensity > 5 x 10 
2 

protons/cm -sec-sr-MeV (which corresponds approximately to 

-3 2 

peak intensity > 1 x 10 proton/cm -sec-sr-MeV for 40 MeV 
proton events). This relation was first observed qualitatively 
on a comparison between the maximum intensities of 10 MeV 
protons and 0. 5-1.1 MeV electrons for 42 events detected by 
the GSFC cosmic ray experiment on the Earth orbiting 
satellites IMP 4 and IMP 5 from May 1967 to October 1972 



-23- 


(Figure 9). A further comparison using the 0.2-1 MeV 
maximum electron intensity from the Caltech cosmic ray 
experiment on IMP 7 versus the 10 MeV maximum proton 
intensity from the GSFC cosmic ray telescopes on IMP 7 
for events detected between October 1972 and September 1974 
(Figure 10) confirms the previous result and gives an 
absolute value for the constant between the intensity of 
0.2-1 MeV electrons (mean energy = 0.4 MeV) and 10 MeV 
protons : 


J Max (0.4 MeV) 

— *\j 100 

J Max (10 MeV) - 
P 


(19) 


Assuming a spectral index of 3 for electrons above 0.2 MeV, 
the ratio at constant energy can be represented by: 


-y ^ !0 2 E °‘ 3 (E in MeV) 

P 

However, this is an average value and one should realize that 
there are fluctuations in the spectrum. 


B. MAXIMUM PROBABLE ELECTRON FLUX AND FLUENCE 

The range of electron energies of interest is _> 500 keV. When the 
-2.7 -3.0 

spectral shapes E * for protons and E ’ for electrons are combined 
with Equation (19) above, we find 


J (0.5 MeV) 
e 


. . . .3.0 /onX 2.7 

M.) 20 , 

\0.5j \10 / 


J (20 MeV) 
P 


333 J (20 MeV) 
P 


Integration over energy gives 


I (> 0.5 
e - 


MeV) 


oo 

-f 


0.5 


J (0 
e 


5 MeV) (E/0.5) 3 dE 


and 


J (0.5 MeV) / 4 
e 



J (20 MeV) (E/20) 
P 


- 2 . 


7 


20 


I (> 20 MeV) 
P _ 


dE 


12 J (20 MeV) 
P 
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Thus , 


I (> 0.5 MeV) 
e — 


7— Vo 1 (> 20 MeV > = 
4 x 12 p — 


7 I (> 20 MeV) 
P - 


( 20 ) 


By analogy with Equation (7) in Section IIIB, we find that the maximum 

probability of encountering a given flux I of > 500 keV electrons at 

4 R is 
s 


P 

max 


<V 


2 x 10 5 
I 

e 


where I is in units of electrons/cm -s. Again, as an example, there is 

6 7-2-1 

a maximum of a 1% chance of encountering an electron flux of 2 x 10 cm s , 

12 

or a maximum probability of 1% of encountering a total fluence of 10 

2 

energetic electrons/ cm during the 16 hour perihelion passage. 


C. MOST LIKELY ELECTRON FLUXES AND FLUENCES 

The proton probability calculations in Section IIIC can be carried 
over directly for electrons if we use the result of Equation (20) and 
let N o (e's >_ 0.5 MeV) = 7 N (p's >. 20 MeV). Thus Figures 7 and 8 can be 
applied to energetic electrons by multiplying the flux and fluence scales 
by a factor of 7. 
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V. SOLAR NEUTRON FLUXES AND FLUENCES 

The estimation of the solar flare neutron flux and fluence near 
the Sun is difficult due to the fact that a neutron flux of solar origin 
has not yet been unambiguously detected. This is understandable, since 
all experimental attempts to identify a solar neutron flux have been 
made in Earth satellites or balloons, where two physical constraints 
make detection of solar neutrons an extremely difficult problem. First, 
the approximately 12-minute half-life of neutrons reduces greatly the 
flux of low energy neutrons which could survive to 1 AU. Secondly, 
local neutron background generated in the spacecraft — mainly the cosmic 
radiation — have placed severe limits on the sensitivity for detection. 
Indeed, only upper limits on neutron fluxes have been reported (Lockwood 
et al. , 1973; Kirsch, 1973). Figure 11 illustrates the effect of the 
decay of neutrons of selected energies as a function of radial distance. 
The measurement of this dependence has been proposed as a means to 
identify solar neutrons and to separate the solar from background neutron 
fluxes (Anglin et al., 1972). Two components of solar neutron fluxes are 
expected to be present, namely: 

A. Quasi-steady solar neutron production which, though primarily 
from extended areas of solar activity, may also exist at lower 
levels and energies over the local photosphere (mainly 10 keV 
to 1 MeV neutrons) , and 

B. Impulsive neutron production from solar flares associated with 
collisions between accelerated charged particles and the 
chromospheric and coronal material (mainly 1-100 MeV neutrons). 

Because neutrons have long mean free paths for propagation through 
spacecraft materials and can produce interstitial damage in semiconductor 
devices, we must decide whether neutrons could represent a significant 
hazard for the proposed Solar Probe missions. In this Section, we 
estimate the neutron flux and fluence possible from solar flare events, 
which undoubtedly have higher intensities and energies than do neutrons 
from quasi-steady production mechanisms (see A above). 
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Two approaches to estimating neutron fluxes from flares are 
possible: 

1) The upper limit for production of secondary nuclei produced 

12 3 3 

in flares could be used for large flares (H , H , H , He , etc.), 
or 

2) evidence from nuclear gamma rays could be invoked. In this note 
we follow the latter approach which was worked out in principle 
for the large flare of August 4, 1972, by Ramaty, Kozlovsky and 
Lingenfelter (1975; referred to as RKL). They calculated both: 

a) the neutron production spectrum, and 

b) the yield of gamma rays due to neutron capture on protons 
(2.2 MeV gamma ray line) and excited states of nuclei in 
the solar atmosphere (e.g., the 4.4 MeV carbon line), 

assuming unit values for the "emission measure" for the flare, 

« that is, the product of the density in the interaction region 

and the intensity of high energy primary protons. The observa- 
tions by Chupp, Forrest, and Suri (1972) of the gamma ray flux 
provide the normalization factors necessary to calculate the 
expected neutron flux for the August 4, 1972, solar flare. 

For the calculation of the neutron flux as a function of distance 
from the Sun, a neutron emission flux is used which is of the form: 

E n < 30 MeV 

-2 

for En > 30 MeV 

-2 

i.e., independent of energy up to 30 MeV and then falling as E (after 
Figure 7 of RKL) . 

For the 4 August 1972 solar flare the normalization factor is: 

25 

= q (E=0) • ftgNp ( > 30 MeV) • = 6 x 10 neutrons/MeV sec 



q (E=0) is the neutron production rate for unit emission measure, 

—18 —I -1 

taken to be 2 x 10 MeV sec , (after Figure 7 of RKL) 


where 
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n N (> 30 MeV) is the emission measure, calculated by RKL on the 
H P 43 

basis of the gamma-ray observations to be ^5 x 10 for the 

solar flare of August 4, 1972, and 

P is the probability of neutron escape from the production region, 

E 

taken to be 0.6 independent of neutron energy (Wang and Ramaty, 
1974) . 

These results were taken from the RKL model calculation for the thin target 
model with a power law primary spectrum but the results are not changed 
significantly (factor of ^5) for other models calculated by RKL. 

The solar neutron flux as a function of neutron energy and distance 
from the Sun, is then calculated from 

J (E n ,R) = q (E n ) * Ps (E n ,R) • R ~ 2 


where 


Ps (E n ,R) is the probability that a neutron of energy E n will 

survive to travel a distance R. This function, divided by 
2 

R is plotted in Figure 11 for neutron energies of 10 keV, 

1 MeV, and 100 MeV. 

In Figure 12 is plotted our calculation of the expected neutron energy 

spectrum for the August 4, 1972, flare at distances from the flare site 

of 5, 10, and 20 solar radii and 0.2, 0.5, and 1.0 AU. In Figures 13 

and 14, the integral neutron flux (integrated over neutron energy to 
3 

10 MeV) is plotted for R = 0.01 to 1.0 AU (Figure 13) and R = 2 to 
50 solar radii (Figure 14). (The major contributions to this integral 
are from neutrons of energy 1 MeV to 100 MeV. ) 

The extrapolation of these results to flares other than the 4 August 
1972 flare is, at best, tenuous. With neither neutron nor gamma ray 
observations, the only other normalization factor available is the high 
energy proton flux. 

Of 25 solar flares defined as "major" by King (1974) which occurred 
during solar cycle 20 from 1966 to 1972, the August 4, 1972, flare was 
the largest. Total fluences of protons with energy above 10 MeV ranged 
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from ^10% to 0.1% of the August flare fluence. Presumably the expected 
neutron flux may be scaled from flare to flare in proportion to the 
observed proton fluence. Thus one may naively expect an average "major" 
flare to contribute a neutron flux of a few percent of the values 
calculated here. Finally, due to the fact that the neutron pulse from 
a flare is so short (on the order of an hour) any observer will either 
see the flare and be subject to the full neutron fluence or will be on 
the wrong side of the Sun and see nothing. Assuming an isotropic neutron 
angular emission distribution, any observer will see 'v half of all neutron 
events. 

In summary, for the spacecraft on a trajectory which passes within 
10 R g of the Sun, the neutron fluence for a flare on the "visible" side 

of the Sun of the magnitude of the August 4, 1972, event, and with 

3 6 7 2 

neutron emission occurring for VL0 J sec, will be VL0 D - 10 X neutrons/cm 

in the energy range 1-100 MeV. 
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VI. CONCLUSIONS 

In the present concept of the Solar Probe mission, an eccentric, 
low perihelion orbit is attained by a Jupiter gravity-assist maneuver. 
The near- Jupiter trajectory is similar to that for the Solar Polar 
mission for which the predicted energetic particle populations are 
given in Table 3. For ease of the implementation of the Solar Probe 
mission, we require that there be no more than a 1% chance that the 
peak flux and fluence of energetic particles near the Sun exceed that 
of the Jupiter flyby. 

TABLE 3. PREDICTED ENERGETIC PARTICLE ENVIRONMENT FOR 
JUPITER PHASE OF THE SOLAR POLAR MISSION 
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2 

Note from Table 3 that at Jupiter I (> 0.5 MeV)/I (> 20 MeV) = 10 , 

e p 

whereas the solar particle energy spectra yield a ratio of only 7 (see 
Equation (20)). Thus, energetic electrons should not be a problem 
either. 

For the Solar Polar mission, the predicted flux of 1 - 3 MeV 

-2 -1 

neutrons from the radioisotope thermal generator is 80 cm s . From 

Figure 12, we see that for the very large flare of 4 August, 1972, the 

-2 -1 

neutron flux in this energy range was perhaps 1400 cm s at 5 R^. 

Thus, in this energy range, we expect the solar neutron flux to exceed 
the background level from the power supply only a few hours a year. Thus, 
although it is unlikely that solar neutrons would be hazardous to the 
health of the spacecraft, there is a severe problem in limiting the 
neutron flux from the power supply to levels which allow solar neutrons 
to be measured. 

One might conclude that flying by the Sun at 4 is less hazardous 

than flying by Jupiter. A solar flyby is not without risks, however. 

First, we could be unlucky and meet an extremely rare, but large flare. 

36 

About one flare per ten or twenty years may emit on the order of 10 

energetic protons, which, according to Equations (15) and (16), would 

7 -2 -1 

yield fluxes and fluences at the Solar Probe of 6 x 10 cm s and 
12 -2 

10 cm , respectively. Such a flux would far exceed the peak flux of 
the Jupiter encounter; the fluence would be approximately the same. 

Second, some of the assumptions made in Sections III and IV may not be 
correct. For example, we assumed that particles travel from the Sun to 
1 AU with no change of energy. If there were some mechanism which 
causes an adiabatic energy loss near the Sun, the flux and fluence at 
a given energy could be very much higher near the Sun than we have 
calculated. Third, although low energy particles (v 1 MeV) are not a 
severe hazard for the internal components of the spacecraft, they may 
be a serious concern for surface elements such as thermal control 
surfaces and sensor surfaces (optical, etc.). An important difference 
between low and high energy particles from solar events is that at low 
energies the events may last for many days with very long injection 
times. Therefore, the probable flux of low energy particles may not 
simply scale as the observed energy spectrum near the peaks of events. 
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Finally, although, with proper shielding, the Jupiter encounter radiation 
levels may not be physically damaging to the equipment, they may very 
well be a strong source of interference. For example, special precautions 
may be necessary to prevent energetic particles from causing unacceptable 
background counting rates or currents in the plasma and imaging instruments 

In summary, although we must continue to improve our models over the 
next few years, a first quick look at the problem indicates that energetic 
particles are probably not a safety hazard for a 4 R g flyby of the Sun. 
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Figure 3. Measured spectral index y as a function of heliolongitude. 
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Figure 4. Time histories of the number of proton events, the number of 
neutron events, and the sunspot number. 
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Figure 5. Sunspot number plotted versus proton event frequency 
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Figure 6. The relation of the number of events per unit intensity as a 
function of intensity for data from 1967 to 1972. 
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Figure 7. The probabilities of the spacecraft encountering fluxes of 
given levels for sunspot numbers of 100 and 200 per year 
and for immediate escape and bubble transients. The 
steepest line is the maximum probability as defined in 
Section IIIB. 
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MAX. INTENSITY, 10 MeV PROTONS/CM - SEC-SR - MeV 

Figure 9. Maximum intensities for electron and proton events are plotted to show 
constant ratio. These data are from the GSFC group from 1967 to 1972. 
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Distance from Sun (A.U.) 

Figure 11. The dependence of the normalized neutron flux 
as a function of radial distance from the Sun 
for three decades of neutron energies. 
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Figure 13. The neutron flux density integrated over neutron energy to 
1000 MeV versus solar distance. 
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